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Abstract 

The decay of highly excited massive string states in compactified heterotic string theories 
is discussed. We calculate the decay rate and spectrum of states carrying momentum and 
winding in the compactified direction. The longest lived states in the spectrum are near 
BPS states whose decay is dominated by a single decay channel of massless radiation which 
brings the state closer to being BPS. 



1. Introduction 

In general, highly excited strings are unstable and will decay. The decay process poses 
interesting questions. What is the lifetime of the unstable string? What is the decay 
spectrum? Is the decay dominated by massless radiation or the splitting of the string into 
two massive strings ? These questions have been investigated for various string theories in 
the past. 

The decay of highly excited open string states was investigated in Q] @ @ [|J • The lifetime 
of a string state of mass mn was determined to be T ~ -4- — . This is consistent with the 

o u gjm 

facts, that an open string has a constant splitting probability per unit length and the size 
of an excited string grows linearly with the mass. 

The decay of closed strings is expected to be qualitatively and quantitatively different. 
Although a generic closed string state is expected to behave similarly to the open string 
H||, the fact, that a closed string can only split (in the absence of D-branes) when two 
distinct points of the closed string are coincident, leads to interesting new phenomena. In 



a series of papers [f?[|]q]P| [Pl||p.H , Iengo and Russo found particular excited string states, 
whose lifetime grows with mass like T = ^n^o- Therefore, for large masses, these states 
are very long lived. Furthermore, the decay into two massive strings is exponentially 
suppressed. The authors explained their results, using a semiclassical argument, relating 
the long lived quantum string state to a classical solution, representing a rotating ring. 
The string never self-intersects in this configuration, so that, the dominant decay mode is 
through the emission of massless radiation. 

Supersymmetry provides a completely different mechanism for the stability of massive 
string states. A BPS state is annihilated by some supersymmetry charges and the super- 
symmetry algebra implies a relation between the mass and the charge of the state. Such 
BPS states are generically stable against decay (they can decay into other BPS states 
at special points of marginal stability). We ask, whether a state being near BPS, in the 
appropriate sense, would also be long lived. 
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The specific string theory, we are considering, is the heterotic string compactified on a 
circle of radius R. Due to the difficulty of constructing the vertex operators for general 
massive string states, we use a very specialized state, involving only level 1/2 and level 
1 oscillators (i.e. states on the leading Regge trajectories). The set of states is uniquely 
parameterized by the mass and the left and right charges. We compute the decay rate and 
the spectrum for such states. 

In particular, the compactification allows us to study massive states which are close to 
the BPS bound, since a state with many left oscillators, but few right oscillators can be 
constructed. We show, that these states are, indeed, made long lived by supersymmetry 
and that, the dominant decay mode is via the emission of massless radiation. In addition 
to the BPS states, satisfying rriQ = k\, the heterotic string also has extremal states, which 
satisfy ttiq = k\. We contrast the decays of 'near extremal' to the 'near BPS' states. 

Let us briefly review several alternative methods that can be used to analyze the decays 
of massive string states before moving on to the approach used in this paper. 

First, the three point function (i \ V(0) \ f) can be used to directly calculate the 
amplitude for a specific decay, where the string state i changes to /, by the emission of 
a second string state, represented by the vertex operator V. This method is particularly 
efficient for calculating decays involving massless radiation and extracting the spectrum of 
the massless radiation. 



A variation of this method considers averaging over initial states [12 1. The calculation 



of |I^1|f?lPl starts with the three point function and considers massless emission only, 
but sums over all final states with a fixed mass and averages over the initial states with 
a fixed mass. If the masses of the initial string state and the massive decay product are 
large, the resulting decay amplitude can be evaluated using a saddle point approximation. 
The resulting massless emission spectrum is that, of a black body, with temperature and 
greybody factors depending on the particular string which is considered. 



A second method ||15|| ||16|| extracts the mass shift for a massive string state, from the 
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residue of the double pole in the s-channel of the one loop scattering amplitude of four 
massless states. The analytic continuation in external momenta is well defined for the four 
point function and the double pole can be unambiguously identified. 

In this paper, we opt to use the optical theorem, which relates the imaginary part of the 
one loop two point function, i.e. mass shift Sitiq of a string state of mass mo, to the total 
decay rate V of this state 000. 



r = Im{8ml) 

The mass shift Sm 2 is given by the one loop string amplitude with two vertex operators 
associated with the initial massive string state inserted. 

5ml = M = g 2 s j d 2 r(v(0) J d 2 zV\z,z)) (1.2) 

where g s is the string coupling constant. It is possible to extract partial rates for the 
decay into two string states of a given mass. However, the rates are always inclusive, 
in the sense that one automatically sums over all decay products and polarizations of a 
given mass. The final computation is performed by expanding the amplitude in powers 
of e 27 ™ T , e lz ,e _27riT , e~ lz \ picking out the contributing monomials using mass and charge 
conservation; analytically evaluating the integrals; and numerically evaluating monomial 
coefficients. 

The relation of black holes and excited fundamental strings goes back to [[19] []20f| ||21|1 [p2 



The exact identification of the (small) black hole entropy and BPS string states 
23(1 [f2~lH [p5| [ p6| Wh was achieved using a' corrections to the supergravity action |28| |29| [j30 |. 



BPS black holes have zero temperature and are stable, whereas non BPS black holes do 
decay via Hawking radiation. While it is interesting to pursue the relation of the decay 
of perturbative string states and Hawking radiation in light of the string/black hole corre- 
spondence, we will limit ourselves to the perturbative string in this paper. For some recent 



work on the relation of absorption in perturbative strings and black holes, see [31 



The organization of this paper is as follows. In section 2, we discuss the spectrum of the 
heterotic string, the condition for states to be BPS, or near BPS, and construct the vertex 
operators of the massive string states that we will use in the rest of the paper. In section 3, 
we calculate the two point function for these states and extract the decay rate. In section 
4, we present the results for the decay of near BPS states. In section 5, we address the 
decays of other string states, which are not close to BPS. We close with a brief discussion 
of our results in section 6. 



2. BPS and near BPS states in heterotic string theory 

We consider the compactification of the heterotic £0(32) string on a circle of radius R. 
The S 1 direction is denoted by X 9 and the non compact spacetime directions are given by 
with fx = 0, 1, • • • , 8. The leftmoving modes make up the ten dimensional superstring 
and the rightmoving modes make up the 26 dimensional bosonic strings. The sixteen extra 
chiral bosons are compactified on a SO (32) lattice. However, we choose not to turn on 
any Wilson lines along the S 1 and to have no oscillators along either the S 1 or heterotic 
lattice directions. 

Since the right and left momenta, on the circle, break up into momentum proper, /, and 
winding, m, the U(l) charges are of the following form (in the absence of Wilson lines)0, 

/ mR , / mR 

The physical state condition for the left and the right movers 

1 / TflR 1 

L | state) = (-(k^k^ + (— H ) 2 ) + N) | state) = - | state) 

2 R 2 2 

L | state) = (-(/^/c M + (4 - ^) 2 ) + N) | state) = 1 | state) 
2 R 2 



(2.2) 



lead to the following expressions for the masses 

ml = k\ + 2N- 1, ml = k 2 R + 2{N - 1) (2.3) 



1 In these formulas and in the rest of the paper we have set a' = 2. 



From which the level matching constraint follows 



lm + N-N+^ = (2.4) 

The supersymmetry algebra of the heterotic string with one compactified direction is 

{Q a , Ql) = 2(Cr^V + ) ab k» + 2(CT 9 V+) ab k L (2.5) 

Where k^ is given by ( |2.1j ). The supersymmetry algebra ( |2.5|) implies that massive states 
can preserve eight of the sixteen supersymmetries if they saturate the BPS bound 

ml = k\ (2.6) 

It follows from ( |2.3|) that a perturbative BPS state has N = h. For notational convenience 
it is useful to introduce a shifted leftmoving excitation number K defined by 

K = N — - (2.7) 
2 v ; 

A BPS state hence has K = and non-BPS states have K > 0. It follows from ( J2.4[ ) that 
if momentum and winding are nonzero a BPS state can have a large N, corresponding to a 
highly excited string. Since BPS states are exactly stable, in order to discuss the lifetime 
and decays of states, we have to move away from the BPS bound. In the following we 
define a measure of a state's closeness to BPS. 

2K 2(N-1) , \ 

nx = — , n 2 = -i g-^ 2.8 

rrig rriQ 

A state is close to BPS if n\ — > (which implies from ( |2.3| ), that there are few left 
oscillators). For fixed m, the more excited a string is, the larger n 2 (since from ( |2.3| ) this 
implies the largest allowed number of right oscillators). For example, a state with large /, 
small K, and zero winding, m = 0, will have n\ ~ and is therefore near BPS. However, 
in this case, k^ = kn and so, n 2 ~ 0, meaning iV is small and the string is not excited. 



n 2 



highly excited 
rightmovers 



BPS 



extremal 



highly excited 
leftmovers 



n 1 



Fig. 1: Range of BPS and extremality measures n\ and n% 
On the other hand, / = m with large / and small K would yield a state both close to the 
BPS bound and highly excited. It is this last type of state which most clearly exhibits the 
anticipated long lifetime. 

A state which saturates the BPS bound (|2.6| ) has equal mass and charge. In the circle 
compactification there is a second U(l) charge kn and hence, there is another limit in 
which mass is equated with this charge, the so called 'extremal' limit 



2 _ j, 2 
m _ K R 



(2.9) 



The condition for a highly excited extremal state is that — > and ri\ — > 1. This has 
the same form as the BPS condition with the role of n\ and ni exchanged. Note that this 
bound is, however, unrelated to supersymmetry. 

2.1. States and vertex operators 

In order to calculate the decay rate and lifetimes one has to choose a particular state and 
compute the related vertex operator. In the lightcone gauge, the massive spectrum is easy 
to construct |32| |33[] ; solving the covariant physical state conditions for general massive 
states is quite complicated. We consider special states and polarizations to simplify the 
calculation. It would be interesting to generalize the calculations in this paper to more 
generic states. 

6 



The left moving part of the state is taken to be 

^u 1 ...u K b\a'L\a'L\---a u _ I l | 0) (2.10) 

2 

Physical state constraints of Virasoro currents and supersymmetric currents L\ \ states) = 
and G1/2 \ state) = necessitate that ^ Vx ... VK is transverse to the momentum in all 
indices. G1/2 | state) = also insures symmetry under any 2 index exchange. Finally, 
L2 I state) = and G3/2 | state) = guarantee tracelessness in any pair of indices. 

r^ 2 ...„ K = o, ^w.., K =o (2.11) 

The polarization tensor is normalized as ^ lJLVl --- VK ^ Ul '" UK = 1. 
The corresponding vertex operator in the —1 picture is given by 

V[~ 1] = -^^...^e-^dX^dX" 2 ■ ■ ■ dX VK (2.12) 



In a one loop calculation the vertex operator in the zero picture instead of (|2.12|) is needed. 
The picture changing operator is given by 

E(z) = e^dX^(z) (2.13) 

and the zero picture operator is given by 

V(°\w) = lim E(z)V^(w) \ {z _ w) o (2.14) 

i.e. one picks the term proportional to (z — w)° in the OPE of the vertex operator in the 
—1 picture and the picture changing operator. The normalized vertex operator is 
then given by 



v ( 0) = v K ^ Ul ,„ VK ^dip Ul dX V2 dX V2 ■ ■ ■ dX UK + ■■■ (2.15) 

y/{K-l)\ 

where the dots denote terms which vanish in the two point function because of the sum 
over the spin structures. 



The rightmoving part of the state is taken to be 



U... Pi ,a p l 1 ---a p _\ \k) (2.16) 

and the associated normalized vertex operator is given by 

Vr = -^U... Pfj BX^ ■ ■ ■ dX"* (2.17) 

This state satisfies the physical state conditions if the polarization tensor £ is transverse, 
symmetric and traceless. The polarization tensor is normalized as £v 1 —vf } £ Vl "' l/fi = 1- 
Note that we, furthermore, simplify the calculation by imposing the conditions that the 
polarization tensors £ and £ are orthogonal and do not have any leg in the S\ direction. 

^i-^K+iC 1 ^...^ = 0) &-1/K-+1 = 0, C9-^ = (2-18) 

The complete vertex operator is then given by 

V(Z,S,k,k L ,k R ) = v[°\z)V R (z)e ik ^+ ik ^+ ikRXR (2.19) 



3. Calculation of the two point function 

The optical theorem relates the decay rate to the imaginary part of the mass shift at one 
loop. The mass shift is given by the one loop string amplitude with two vertex operators 
( p.!9|) insertions. 

Sm 2 =A 2 = g 2 s J d 2 r(v(0) J d 2 zV\z,z)^ (3.1) 

where the z integral is over the usual fundamental torus domain, i.e. the parallelogram 
stretched onto the unit vector along the real axis and the r = n + ir 2 vector, and the r 
integral is over the fundamental domain of SL(2, Z), the conformal transformations of the 
torus. 
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In order to calculate A 2 we use the vertex operator (|2.19|) and evaluate the various 
contractions. 

A 2 = c'g 2 s K(£-0(t-0 J d 2 TT 2 J d 2 z(e lk » x " (0)e- lk » x " (z,z)}(e lkLXL (0)e- lkLXL (z)} 
(e lkRXR (0)e- lkRXR (z)}(dX(0)dX(z)} K - 1 (BX(0)dX(z)f 

Y,^Z v ^d^)^(z)) v Z het (f)Z, 

(3.2) 

where we used the fact that the vertex operators are normalized and the orthogonality 
of the polarization tensor to eliminate right with left contractions. The constant d is an 
overall normalization which is independent of the initial state. Since we are interested only 
in the relative comparison of rates and lifetimes we will set d g 2 = 1 when evaluating the 
amplitudes. Z v is the partition function for fermions with spin structure v\ Zh et is the 
partition function for the bosons compactified on the SO (32) lattice; Z is the partition 
function for bosons along the first ten dimensions and for the ghosts. The factor of r 2 
is from gauge fixing the conformal symmetry of the torus by choosing the first insertion 
point to be at 0. 

The contractions and contributions to the partition functions are evaluated in Appendix 
A. The fermionic contractions and the sum over spin structures is performed in (A. 3), the 
SO(32) partition function is given in (A. 7), the bosonic partition function is given in (A. 8) 
and the bosonic correlators are given in (A. 11) and (A. 14-15). Putting everything together 
the result is 

n,w 1 2 i=l 

,2 





r) 


01(0 | 


r) 



2 1 \ I / 

(3.3) 

where z = z\ + iz 2 . The sum over the integers n, w represents the summation over the 
discrete loop momenta = n/R + wR/2 and pu = n/R — wR/2 in the compactified 
direction, which also makes them the left and right charges of one of the decay products. 
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The constant d is a normalization which is independent of all the parameters. The plan 
to evaluate the amplitude is as follows. First, we expand the amplitude in infinite series 
in powers of q and e lz and their complex conjugates. Secondly, we perform an integration 
over the n and Z\ which yield constraints on the exponents of q = e 2 " E%T and e lz and 
their complex conjugates. Next, we cut off the infinite series to obtain finite polynomials, 
by relating the exponents to the masses of the original state and its decay products, and 
applying conservation of mass and charge. More precisely, not only is the series cut off, but 
a finite list of the exact powers of q and e lz and their complex conjugates which contribute 
to the amplitude is compiled. The integrals for all the monomials, over r 2 and Z2 are, 
then, performed analytically. Unfortunately, the final computation of coefficients of all the 
monomials can only be performed, for specified initial states, using Mathematica for the 
power series expansions of the expressions. 

Using the formulae given in appendix B the holomorphic part of ( |3.3|) can be expanded 
as a power series in q and e iz 





r) 


0[(0 


r) 



2K 



K ~ X 1 fK-l\ rK G {:U) 



* a (e lz ) b 



r=0 v ' a,b 

The antiholomorphic part of the amplitude (|3.3| ) can be expanded as 
Using (|3.4| ) and (|3.5| ) the two point amplitude A2 becomes 

K-l N -(r+s+|) 



(3.5) 



a 2 = c' 9 2 s ky: d*r /rf^EE 



r 2 v '"' 2/ f K-l\ f N 



(4 7r )r+s I r I \ s 
r=0 s=0 v ' v 7 v 



^ C K,r C^' s q a+ ^ p2L q c +y 2 Re iz{ - b ~ kLPL "> e~ iz{ - d+kRPR) 



(3.6) 
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The integration over the modulus r in (|3.6| ) is over the fundamental domain of SL(2, Z). 
However, as explained in the next subsection or in [BJ , if one is interested in the imaginary 
part of A2 only (as we are), the fundamental domain can be replaced by a regular strip 
—1/2 < n < 1/2 along the entire length. 

The integration over n and Z\ then take the form: 

• 1/2 r 2-w 
1/2 



/1/2 />2tt 
d Tl e 2niT ^ nw+a - c \ / cfe ie * 21 
-1/2 JO 



(— nm — wl-\-b— d) 



(3.7) 



so a non-zero amplitude requires the a — c + nw = and b — d — nm — wl = 0. After 
performing the integrations over ui and 7*1 the two point amplitude becomes 

K-l N 



a 2 = wrfjr / dr 2 f dV2 y, E E ( A ' r " 1 ) ( s , ) E c S r cf + 

' ^ n.iu r=0 s=0 ^ / \ / „ p. 



s 

nw ,b—nm—wl 



_-(r+»+§) „ , , ™; , 

_[_2 e -27rr 2 (2a+(n/i?+'u;i?/2) 2 ) e -i/2(2&-2nZ/fi 2 -i(;mi? 2 /2 nm wl) ^ 27TT2 ^2 

(47r)r+s 



(3.J 



S.i. Calculation of decay rate 
After a change of variables 2tct2 = t and 1/2 = ty the integrals in (J3.8|) are all of the 



following form 

&2 



[ dt [ dyt <*-^ e -t{<y 2 -y{<-rn\+™i)+™i) ( 3 9 ) 

Jo Jo 



where a = — r — s — 5/2 and is negative. 

The choice of variable names mi, m2 is not accidental. By considering a Schwinger 
parametrization of a Feynman integral for a loop amplitude, we can convince ourselves 
that the coefficients of powers of y are associated with the masses of the decay products 
as above 0,0,0. For completeness, the argument is reproduced in Appendix C. 

A comparison of (|3.9|) with (|3.8|) gives the following expressions for the masses mo, mi, m2 



'l — n (m — w)R\ 2 n , , _ to -\rv\ 

n>- = \— + - + 2 ( a + b + K ) ( 3 - 10 ^ 

, n wR\ 2 
mi=l R + —) +2a 

11 



Note that a comparison with the mass-shell condition (2^3) leads to an identification of 
the momentum ki t L, k 2 ^ as well as the oscillator level Ki, K 2 of the two decay products 
with mass m± and m 2 respectively. 

fl — n (m — w)R\ 
ki,L = k L -PL=(— 5 - + - tt^— , K 1 = a + b + K 

p 27 (3-11) 

k 2 ,L=PL= [^ + —), K 2 = a 

where k^ is the leftmoving compact momentum of the initial state and pl is the compact 
loop momentum. As we shall see later the power series expansion of the amplitude limits 
the range of b to be —a — K < b < —a. 

The optical theorem relates the imaginary part of the loop amplitude to the summed 
squares of the tree three point amplitudes. The loop amplitude is formally real but the 
imaginary part comes from the analytic continuation of the amplitude into a region where 
it is divergent. From ( |3.9|) we can see that the integral is divergent when the polynomial 
coefficient of t in the exponent, 

p {y) = m W ~ y(ml ~ mf + mj) + m\ (3.12) 

is negative, which occurs between the roots of this polynomial, 



= < ~ ™\+ r4 ± J_ Uml-ml + mlY _ 
y 2m 2 m o]l 4mg 2 v ; 

This restricts the domain of integration in the y variable. Note that the roots are real 
and the imaginary part of A 2 exists if and only if mo > m\ + m 2 . 

Let us now turn to the analytic continuation of the t integral. The two divergences to 
deal with are a pole at t = and an essential singularity at t = oo. We will analytically 
continue in P and a. The integral is well defined for Re[P] > 0, i?e[a] > 1. We make the 
substitution u = tP. The t integral is then: 

exp- u u a - 1 du= — I ex^- u u a - l du- — I exp- u u a_1 dw (3.14) 

12 



pa / 1 pa / 1 p 

r 'Pe r JO r 



Pe is the variable lower limit inherited from the r domain. In a moment we will see that 
the exact lower limit does not matter in the calculation, because the second integral is 
real for all values of P and a, and therefore, the complicated r domain can be replaced 
by a strip extending down to 0. (It may worry the reader that this integral domain 
now includes infinite copies of the fundamental domain and should diverge by modular 
invariance. However, modular invariance has, already, been broken in our calculation 
because we do not keep the entire modular invariant expression for the amplitude, but 
instead pick out only a finite number of monomials.) More explicitly, if we analytically 
continue the first integral to negative P we obtain: 



, exp u u a du = - — — — T(a) = - — — — — r : (3.15) 

-l) a (-P) a J P (-P) a K 1 {-P) a T(l - a) sin(aTr) v ! 

where we used that a is a half-integer (see (|3.8|) ).i The sign in (|3.15|) is chosen to give a 



positive mass shift. Now the expression can be analytically continued for negative P and 
negative a. For such values the expression is pure imaginary. Now let's see that there's no 
imaginary contribution from the second integral in ( |3.14]) . Let P take its negative value 
and make the variable change : u — > —u (or equivalently make the u = tP substitution on 
the — e integral before doing any analytic continuation in P) 



( _i)«-i_L_ f P \^u ua -i du = __t_ f P£ y_^_ = e <*Y [- p £) n 

{ } (-P) a Jo (-P) a Jo ^n\u^ Z-^ n \(n-a) 

(3.16) 

This expression is convergent with no poles and it is still so if analytically continued 
for negative a. Since this integral has no imaginary part and the first integral had no 
dependence on the lower bound, the fundamental domain of r can be replaced by a strip 

2 Note that, were the number of compactified dimensions even, a would be an integer and the 
analytic continuation would need slight modification. 
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extending all the way down to zero. All that's left is the y integration. 



J (-P)- a dy = J m-^(y - y_)-«(y+ - y)-«m^(y + - y,) 1 " 2 _ ^ (3.17) 

where the roots, y+ and y_, were defined in (|3.13| ) and a = — r — s — 5/2. 
Combining with the t integral we obtain 



, / x . 2 2r+2s+6 ir T r + s + 7/2 )/ (ml - m? + m\ 2 2 \ r+s + 3 

Im(a 2 ) = 6(m - mi - m 2 ) ; ^4 - — ^ — mi ' 

m T(2r + 2s + 7) V 4m 2 , 



(3.18) 

Where the 0-function enforces the on-shell condition mo > mi + m 2 . Applying ( |3.18j ) 
to the terms appearing in the summation ( |3.8| ) gives the imaginary part of the two point 
amplitude 



Im(A 2 ) = c'glK E £ £ £ Q ^ " ™i " ^) ( ~ ) ( , ) ^2 21 /^ 

r=0 s=0 n,u> a, 6 

xC r,K C s,N Tjj+r + s) 1 / (mg-mj + m|) 2 2 f +r+s 

a>6 a+ „ w ,6- n m-»i r ( 7 + 2r + 2s)mo V 4m 2 , 2 / 

(3.19) 

The total decay rate is given by 

r = Im(5m 2 ) = Im(A 2 ) ^ 
2mo 2mo 

The lifetime of the unstable string state is the inverse of the total decay rate 

T=i (3.21) 

Note that the summation over the loop momentum n and winding w and a, b sums over all 
particles with mass mi and m 2 given by (|3.10| ). If the on shell condition is fullfilled and 
^I'fc ^o+ntB b-nm-wi * s nonzero then the the original string state can decay in this channel. 



The decay rate in this channel is given by the imaginary part (|3.19| ) for a fixed n, w, a, b. 
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4. Decays of near BPS states 

Massive BPS states are generically stable against decays and hence one has to move 
away from exact BPS states to investigate the lifetime of excited string states. An excited 
string state represented by the vertex operator ( |2.19|) is characterized by the momentum, I, 



winding number, m, and the left and right moving excitation level K and N, respectively. 
The level matching condition (|2.4|) provides one relation among the four parameters. The 
measure of how close a state is to BPS is therefore not unique and depends on the definition. 
In the following we will consider states which have n\ close to and n% close to 1. This 
can be achieved by choosing K to be close to zero and choosing /, m such that is large 
and kR is small, so that the state has a large N. In the following, we also fix the radius of 
the compactification circle to be R = 1. 

It is expected that the states which are closest to BPS are the longest lived for a given 
mass. It is therefore interesting to determine the dependence of the lifetime on the mass 
mo of the initial state. In the following plot, states with K = 1 and pr = 0, were 
considered with p^ varying such that the mass ranges from rriQ = 98 to = 308. 



4.4 4 .€ 4.8 5.2 5.4 5.6 



Fig. 2: Log/Log plot of lifetime T of K = 1 as a function of mass TOq. The 
masses range from 



A x 2 -fit for the mass dependence of the lifetime gives 



T = e- 8 - n (mo) 6 - 06 (4.1) 
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As we shall show below, the exact behavior of the lifetime of a near BPS state in the limit 
of large masses is T ~ m®. 



In the following table we list the decay rates and the dominant decay channel of near 
BPS states with momentum I = 9 and winding number m = 16 for varying K. The 
lifetimes are calculated using ( |3.20| ) setting g 2 c' = 1 to fix the normalization. 



K 


ml 


m, 


n 2 


lifetime 


dominant decay channel 


per cent 


1 


291 


0.0068 


0.996 


8705.5 


K x = 0, m? = 289, m\ = 


99.9 


2 


293 


0.013 


0.996 


4376.2 


K x = l,mf = 290.7, = 


98.3 


4 


297 


0.0269 


0.996 


2212.5 


K x = 3, ml = 294.8, m\ = 


95.0 


6 


301 


0.039 


0.996 


1492.1 


K x = h,ml = 298.9, m\ = 


92.0 



Table 1: Lifetime and dominant decay channel of near BPS states 
For near BPS states the lifetime increases when n\ gets closer to zero, which is a measure 
of closeness to BPS, hence the closer to BPS a state (of comparable mass) is, the longer 
lived it is. The decay is dominated by a single decay channel, where one of the decay 
products is a massless graviton and the other is a massive string state where K has been 
decreased by one (subleading decay channels decrease K by more than one unit but give 
a small contribution to the total rate). Hence a massive near BPS state will gradually 
become more BPS by massless radiation. 

The mass dependence ( |4.1|) and the dominance of the massless decay channel can be 
understood in more detail by analyzing (|3.19|) . Since the decay amplitude is symmetric 
under interchange of m\ and m2, for a massless decay we can pick the second particle to 
be massless, i.e. m^ = 0. Conservation of energy implies that the first decay product has 
mass m\ = k\ + 2(K — AK), where AK is a positive integer related to the decrease in 
non BPS-ness. For the massless decays the last factor in ( |3.19| ) then becomes 



ml + ml) 2 m 2^ 3+r+s _ ^AK 2 Y +r+s 



4m,Q 
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mf. 



(4.2) 



The inverse powers of rug in ( f4.2|) suppress terms with r, s > 0, unless these inverse 
powers are balanced by combinatorial factors coming from the power series expansion. 
The combinatorial factors on the holomorphic side are small for near BPS decays because 
they are of order K and K is small, but on the antiholomorphic side combinatorial factors 
can be large since they can be of order N, i.e. of order mg for near BPS states, so we can 
set r = to a good approximation, but different s can contribute equally. Now let's see 
what powers of q, q, e lz , e~ lz we're after. It follows from m2 = and ( |3.10|) that a = as 
well as n = w = 0. This means that in the power series expansion (|3.6| ) we are looking for 
the terms with q°. In this case the mass of the first decay product is given by m\ = rriQ + 2b 
and it follows from mo > m2 that we are looking for a negative power of (e lz ) b in ( |3.4|) . 
From Appendix B we can see that if we keep only terms proportional to q° and set r = 
we get 

oo 

(J2le lzl ) K ~\ie-' z/2 (l - e lz )) 2K (4.3) 
i=i 

We can see that only the first negative power of e lz is present so b = —1. Comparing with 
( [3 .11|) one sees, for a = 0, b = AK, which explains why the dominant decay channel has 
AK = 1. On the antiholomorphic side, from (|3.7p , we know that c = and d = —1, so we 
are looking for q°e lz . Keeping only q° terms from (|3.5|) we obtain 

(- + 24) (2 + 960?) (1 + (-6 + 2N)q)(ie lz/2 (l - e~ lz - qe- 2lz )f N - 2 

oo oo (4.4) 

x ( ( E le ~ l ~ Zl ) + ( N ~ s ^ e ~ l ~ Z + e ") ( E le-** 1 )"'" 1 ) 
i=i i=i 

Expanding the above expression shows that for small s the coefficient of (fe lz goes like 
iV s so small s contribute to the computation at the same order of tuq. (For large s the 
combinatorial factors no longer go like iV s .) Therefore we see that, in the limit of large 
initial mass, the decay rate of near BPS states behaves like Y ~ 1/toq in accordance with 
the numerical result (|4.1| ). 
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5. Decays of other states 

The calculation of the decay rate and lifetimes is not limited to near BPS states. It is 
an important question, whether the near BPS states are indeed the longest lived states. 
To this end, we have calculated the lifetimes for all inital states of a given fixed mass, 
ml = 171. For radius R = 1, there are 134 possible initial states. In the following table, 
we list the ten longest lived states among them. 





kR 


K 


ni 


n 2 


lifetime 


decay channel 


-13 


-1 


1 


0.0117 


0.9942 


1750 


massless 


-13 


1 


1 


0.0117 


0.9942 


1683 


massless 


-13 


-3 


1 


0.0117 


0.9474 


1645 


massless 


-13 


3 


1 


0.0117 


0.9474 


1197 


massless 


-13 


5 


1 


0.0117 


0.8538 


375.8 


light 


-13 


-5 


1 


0.0117 


0.8538 


99.75 


light 


-13 


7 


1 


0.0117 


0.7135 


85.95 


light 


-5 


13 


73 


0.8538 


0.0117 


62.24 


light 


-7 


13 


61 


0.7135 


0.0117 


41.95 


light 


-11 


-3 


25 


0.2924 


0.9474 


38.32 


massless 
















7 


+13 


61 


0.7135 


0.0117 


0.0152 


light 


-13 


+13 


1 


0.0117 


0.0117 


0.0111 


massless 



Table 2: Longest lived states of mass m\ = 171 

The longest lived states are indeed the near BPS ones with massless radiation as dom- 
inant decay channel. A near extremal state appears as the 8th longest lived state, the 
lifetime is however 30 times shorter than that of the logest lived near BPS state. The 
dominant 'light' decay channels consist of BPS states with nonzero winding or momentum 
and a mass close to zero. The states with dominant 'light' decay channels exhibit shorter 
lifetimes compared to the states with dominant massless decays. 

For comparison purposes, we have listed the shortest lived state among the 134 states 
of mass TOq = 171. Note that the state is technically near BPS since ri\ is small, however, 
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it is not highly excited, since n<i is small. Hence, it is the highly excited near BPS states 
which are long lived. 

In contrast to the BPS condition, extremality (|2.9| ), is not related to supersymmetry 
and it is an interesting question whether highly excited extremal states are also long lived. 
Already from table 2 it is clear that extremal states are shorter lived than BPS states (for 
a particular mass). We have calculated the lifetimes of near extremal states with N = 2 
and masses ranging from uiq = 38 to rriQ = 258. 



-2 



-2.5 



4.5 5 5.5 

-3 . 5 

Fig. 3: Log/Log plot of lifetime T of N = 2 near extremal states as a 
function of mass tuq. 

A fit of the data produces a lifetime growing as T ~ m^ 04 , however, note that, compared 
to the near BPS states, the lifetimes of the near extremal states are several orders of 
magnitude shorter. 

A general feature of all decays, whether the initial state is near BPS or not, is that all 
decays are dominated (to at least 99% of the amplitude) by channels where one decay 
product (say the first one) has oscillator level K\ = 0. Note that, in general, the dominant 
decay channel may be not massless, but the mass of the state will come from momentum 
and/or winding, not from left oscillators. To understand this, we must again look at 
the expansion of the holomorphic piece (|3.4|) , but we do not assume a is 0. Remember 
that a = K2 and a + b + K = K±, which imply together with mass conservation, that 
— K < a + b < and a < K, where a is the power of q and b is the power of e lz . Also, 
remember that terms with r > are suppressed. One obvious term one can pick out for 
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r = is (qe~ xz ) K ~ 1 e~' lzk . (This is obtained by picking up from (B.2) only terms from 
second sum with I, n = 1 and only terms with n = 1 in (B.3) .) This term has a + b = —K, 
i.e. K\ = 0, and a = K — 1, as we wanted to argue. Now suppose we try to pick a + b less 
negative. Whether you do this by picking n > in one of the 2K factors in (B.3) or by 
choosing n > 1 in one factor of (B.2), the power of q becomes too large and we are forced 
to use r > 0, which produces suppressed decay rates. 



6. Discussion 

In this paper we have analyzed the decay of highly excited string states for compactified 
heterotic strings. We found that the longest lived states are highly excited, near BPS 
states, for which the lifetime grows as T ~ l/g 2 s (m ) 6 for large masses. The decay of 
such states is dominated by massless radiation, which reduces the non-BPS'ness K by 
one. The end result of such a decay (over several steps) is a stable BPS state. For large 
mass mo, the near BPS states is very long lived. Near extremal states and generic states, 
which have both left and right excitations, are much shorter lived, although the lifetime of 
near extremal states grows with mass. It is an interesting question, whether the decay of 
excited strings can be used for exploring the string/black hole correspondence. To achieve 
this, one has to consider more generic string states, instead of states on the leading Regge 
trajectories. Some steps were taken in this direction, for the comparison of absorption of 
massless scalars by perturbative strings and BPS black holes, in a recent paper and 



it would be interesting to generalize this work to decays of near BPS states. We plan to 
address this question in the future. 
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Appendix A. Collection of formulae for one loop string amplitudes 

The partition functions of the fermionic leftmovers with even spin structure v = 2, 3, 4 
are 

9 2 (0 | r) 4 9 3 (0 | r) 4 fl 4 (0 | r) 4 

r7(r) 4 ?](r) 4 ??(t) 4 

The propagator for two world sheet spinors for the three even spin structures, denoted 
by v = 2, 3, 4, is given by 





r)6>i(0 


r) 


2tt 0„(O 


r) 




r) 



W{Z)1I} V {W)) V = l/*" — nf X/, /,-«, ( A - 2 ) 



Correlation functions for the odd spin structure i/ = 1 vanish, unless there are ten 
fermionic zero modes inserted, which is not the case for our amplitude. 

After summing over spin structures, the correlators involving fewer than four fermionic 
fields vanish. Also, fermionic correlators involving four fermions without derivatives vanish. 
The only nonzero fermionic correlator in our calculation of the decay rate is 

Y,^Z u (r){r{z)d z r(zW(w)d w ^\w)) v = (^V P + ^V A )^(t) 8 (A.3) 

V 

where e„ are phases implementing the GSO projection via €2 = —1, £3 = 1, £4 = —1. 
Calculating (|A.3[) and vanquishing the other fermionic correlators is accomplished by 
applying the following Riemann identity to derivatives of (|A.2j ) and ( |A.1| ) . 

1 4 4 

sE^n^kH-n^M (a.4) 



where 



2 



v[ = -^(-Vl +V2 + V3+ V4), v' 2 = l,(+Vl -V2+V3 + V 4 ) 
V 3 = 2 ( Ul + ^2 - ^3 + U4), ^4 =2(^1+^2+^3- V4) 



(A.5) 



and 



6[(0 I r) = 2nri(rf (A.6) 
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The partition function for the bosonic rightmovers in the SO (32) lattice directions is 



4 



=U6 



Zhet ~ 2 ^ V (f)™ (A - 7) 

The partition function (including the zero modes integral) of the nine left and right 
moving bosons in the uncompactified spacetime directions is 



2 1 1 



The correlation function for noncompact worldsheet bosons X^, \i = 0, • • • , 8 is given by 
(X f "(z, z)X u (w, w)) = -^V*" In ( | E(z, w) | 2 +-^—(z - w - z + w) 2 ) (A.9) 

2 V 07TT2 / 

where the prime form E is defined as 

EM=2 *^mw (A ' 10) 

The correlators involving dX^ can be derived from (|A.9| ) and read 

(dX»(z)dX v (w)) =?Tyd 2 2 ln { dl{Z -^T 1 T) ) +?r 8^ 2 



, / 



(dX»(z)8X»(w)) =<^a/ln I f )) +^ °' 



(A.ll) 



2 Z V v 27T ' V 87TT2 

Using (|A.9|) the following correlation function is also calculated: 



n , 

a 



( "Q e ^x (2i> ,- i)> = JJ | ^ ^ r 'k lk] exp ^ _ _ _ z . + _ i)2 j (AJ2) 

In the 9th direction, compactified along the Si, the zero modes integral becomes a 
sum since the momenta are discrete, so the partition function for the boson along the Si 
direction is 

z Sx = E e^ Tp2 -^ p2 1 (A.13) 
For compact bosons we need a 'crural' version of ( [A. 12| ) (derived in P^f), 

n 
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where p is the loop momentum and ki is the left moving incoming momentum. 
For the antiholomorphic correlators one gets 

n 

(U e***** 2 *)) = exp(-ip hzi) II E(z it (A.15) 

where pi, ki are the corresponding right moving momenta. Note that, since there is 
winding, a compact boson has ki 7^ fcj. We can check, that if multiply (|A.15| ) , (|A.14|) , 



and Zs 1 and replace the sums over pi by integrals, we recover the nonchiral correlator times 



the partition function of one uncompactified boson i.e. ( |A.12|) multiplied by 1/2 } — . 

V v(r)r](T) 

The last piece we need is the ghost partition function 

Finally, to account for the gauge fixing of the conformal invariance of the torus i.e. fixing 
the first insertion point to be at 0, we must augment Z g h with a factor of t%. All the 
factors of r 2 combine to yield r 2 ' . 

For notational convenience we define a single partition function combining the bosonic 
and ghost partition function 

Z = Z st Z Sl Z gh (A.17) 



Appendix B. Expansion of theta functions 

In this appendix we gather the formulae for the series expansion of theta functions and 
related objects in q and e lz . The basic expansion formulae for the theta functions are 

e x {v\r)= i J2(-l) n q^ n -h) 2 e iz ( n -^ 

n 



e 2 (v I r) = ^ 5 3(»-3)%.**(»-3) 



1 2 n2 e* zn 



9 3 (v\r) = J2<l 

n 

4 ( v I r) = ^2(-l) n q? n2 e tz " 



e 

(B.l) 
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From these basic formulas the expansion of the terms appearing in appendix A can be 
derived 

oo 

d z nn{e 1 (^-\r))=Y,Ke lz ) l + E ^ ((e 12 )' + (e^)"') (B.2) 

1=1 l>0,n>0 

and 

i(0|r) ^ ' n m>0 (l-? m ) 3 



9 

and similarly for the antiholomorphic terms 



Appendix C. Identification of Masses in Amplitude by Comparison with 
Schwinger Reparametrized Feynman Integral 

Consider the Feynman integral for the one loop amplitude. For simplicity, let's first 
consider the loop particles to be scalars. 

/OO i i 
dd P-9 27 PT9 2 ( ai ) 
.oo p 2 + mi {p - k) 1 + mf 

where p is the loop momentum, k is the incoming momentum, d is the uncompactified 

spacetime dimension (which is 9 for us). We preform a Schwinger reparametrization and 

obtain 

/OO /"OO />oo 

d d p d<y c / /3e -^(P 2 +^) e - l /3((fc-p) 2 +m?) ( C . 2 ) 
-oo JO JO 

Performing the momentum integration we obtain 

p OO /* OO 2 2 d 

- dry I d/? e -^ m 2-^( fe2 + m ?) e ^w( — - — y (C.3) 
Jo Jo K n + i/3 J 

Finally let's change to variables y = t = i-f + i(3, and rotate the t contour to run 
along the real axis. 



OO p 1 

I / 7T \ ./2 2 / 2 i 2 2\ i 2n 

(it I dy (— ) 2 e ( m ° y -yy m o+ m 2— m i)+ m 2i (C.4) 



o ./o 
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where we used — k 2 = m\ and that the jacobian for the variable change is —t. This 
integral looks quite similar to (|3.9|) , except the power of t in (pT9|) is — r — s — (d — l)/2 
whereas here it is just {d — l)/2. The second issue is that there is no reason to assume that 
the decay products are scalars. Let us, then, consider what happens with some arbitrary 
polarization. Note that the precise polarizations allowed by string theory are not arbitrary. 
Furthermore, the states we are working with have even more specialized polarizations. We 
are not trying to establish an exact correspondence between the string theory calculation 
and the equivalent field theory calculation, but only to extract the masses of the decay 
products. To that end, we just note, that the integral must be a linear combination of 
integrals such as the following 

r p2i(pk)c 1 d d P (c.5) 

J-oo P 2 -m\{p- k) 2 - m\ 
The powers of p must be expanded in components so the expression after Schwinger 
reparametrization looks something like 

d d p d 1 d/3 V ITi}^ e -^(P 2 --^) e -z/3((fc- P ) 2 -^) ( C6) 

-oo JO JO 

where Di are some coefficients, which might depend on k%. The thing to notice here is, 
that though this expression is much messier than the one for scalars, the exponentials are 
unchanged, so performing the gaussian integrals we get something like 

W dry dPe-^<-^ k2+m ^e^(—^—y +9 ° (C.7) 
a Jo Jo + ift 

The momentum integral must be performed in components so the sum is complicated; 

D' a are some messy coefficients, which can be functions of y = /3/(7 + 0) and kf, g a are 

integer- valued. Looking at the gaussian formula, one can be convinced, that the form of 

the exponentials remains exactly the same for each term as in the scalar case. Therefore, 

as in the scalar case, we're justified in associating the masses with the coefficients of the y 

powers in the exponent of the exponentials. 
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